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Abstract 

We study the thermoelectric properties of monolayers of M 0 S 2 and other group-VI dichalco¬ 
genides under circularly polarized off-resonant light. Analytical expressions are derived for the 
Berry phase mediated magnetic moment, orbital magnetization, as well as thermal and Nernst 
conductivities. Tuning of the band gap by off-resonant light enhances the spin splitting in both 
the valence and conduction bands and, thus, leads to a dramatic improvement of the spin and 
valley thermoelectric properties. 


and U. Schwingenschlog 


ogn 


'^m.taliirOGQalumni.imperial.ac.uk 
udo.schwingenschlogl@kaust.edu.sa,-1-966(0)544700080 


1 





I. INTRODUCTION 


fl. 


Being the first truly two dimensional material [l[, graphene has attracted remarkable 
attention, both due to its exotic transport behavior and technological applications in various 
helds j^. Still, fundamental problems restrict its applicability, in particular the negligible 
band gap and weak spin orbit coupling (SOC). These limitations could be overcome W 
monolayer M 0 S 2 , which therefore is interesting for next generation nanoelectronics 
M 0 S 2 combines the honeycomb structure of graphene with a large intrinsic direct band gap 
of 2A = 1.66 eV and a large SOC of A = 74 meV, providing mass to the Dirac fermions 
js-lOj. As a consequence, preliminary results indicate potential in valleytronics, because the 


dispersion can be manipulated in a flexible manner for optoelectronic applications 


10 


13|. 


Spin and valley Hall effects have been predicted in an experimentally accessible temperature 
regime jo], the former arising from the strong SOC and the latter from the broken inversion 
symmetry. 

In addition to the electrical and optical properties. Berry phase mediated thermoelec¬ 
tric effects due to a temperature gradient have been proposed for two-dimensional systems 
ji^ . Orbital magnetic moments, orbital ma gnet izations 


conductivities have been addressed in Refs. 


17- 


15| . as well as thermal and Nernst 


19l | and theoretical models for the thermo¬ 


electric transport have been presented for graphene in Ref. 20|] and for topological insulators 


in Ref. 


2l|. Of particular interest is the tuning of the spin and valley thermoelectric prop¬ 


erties of M 0 S 2 and other group-VI dichalcogenides, where a temperature gradient gives rise 
to transverse spin/valley accumulation and spin/valley current. In graphene this is difficult 
to realize due to the negligible band gap and weak SOC. 

In the present work we quantify the Berry phase mediated thermoelectric properties of 
M 0 S 2 and other group-VI dichalcogenides by deriving analytical expressions for the key 
thermoelectric quantities in the presence of circularly polarized off-resonant light. Gap 
opening by such light has been predicted for graphene and for the surface states of topological 
insulators 22|, and has been conhrmed experimentally for the latter 23|. For graphene the 


chiralities for different values of the frequency have been given in Ref. 2^. Moreover, 
opening of a trivial gap has been reported under high-frequency linearly polarized light 25 1. 
Going beyond these hndings, we demonstrate in the following that by off-resonant light large 
spin and valley thermoelectric effects can be achieved. 
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II. MODEL FORMULATION 


Extending the approach of Ref. 
effective Hamiltonian 


by introducing time dependence, we start from the 


= v{r]a^Il^(t) + aylly{t)) + Acr^ - Ar^s^cr^ + XrjsJ (1) 

in the a;|/-plane in the presence of circularly polarized light, where r] = ±1 represents the 
K- and A'-valleys, respectively, A is the mass term that breaks the inversion symmetry, 
CTj;, ay, and are the Pauli matrices, A is the SOC with real spin index s^, and v = ato/h 
is the Fermi velocity of the Dirac fermions (with to being the nearest neighbour hopping 
amplitude and a the lattice constant). We use the gauge in the two-dimensional canonical 
momentum 11 (f) = p + eA{t) with vector potential 

A(t) = (iAsinOt, AcosOt), (2) 


where D is the frequency of the light and A = Eq/XI {Eo is the amplitude of the electric 
held, E(t) = —dA{t)/dt). We have A{t + To) = A{t) for Tq = 27r/D. The plus/minus sign 
in Eq. (2) refers to right/left-handed circular polarization of the light. 

The effect of off-resonant light can be described W a static Floquet Hamiltonian 22| . 


which yields excellent agreement with experiments 


23|. A static approach is satished for 


low intensity [evA <C Kl) and high frequency {to <C Kl) light, which does not directly excite 
electrons but effectively modihes the band structure through virtual photon absorption and 
emission processes. We arrive at the effective Hamiltonian (see the Appendix) 

= viWxPx + (TyPy) (A 7]An)a^ - Xps^a^ XrjsJ, (3) 

where Aq = is an effective energy term representing the circularly polarized 

off-resonant light, which essentially renormalizes the mass of the Dirac fermions. Similar 


approaches have been used for describing gapped systems such as silicene 


26| and disordered 


topological insulators [27[. Diagonalization of the Hamiltonian leads to the eigenvalues 


Ep"^ = sr]X + ( \/ (vHk)'^ -|- (A -|- rjAfi — Xrjs)'^ 


( 4 ) 


and eigenfunctions 


T"’’^(k) = 


/ 


^ik^x+ikyy 

LxLy 


vhke-^^v 


y (vhkP+ |^-A-»7An+77sA+C-^ (D/ifc)2 + (A+An-Ar;s)2’j 
—A—r;An+r;sA+C ('ufife)^ + (A+An-Ar;s)2 
^ {vhkP+ |^-A-»7An+77sA+Cy^ {vfikP + {A+An-Xr]sp'\ j 


( 5 ) 
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Here C = =tl denotes electron/hole bands, s = ±1 stands for spin up/down, and tp = 
taii~^{ky/kx) with kx = k cos ip, ky = ksinp, and k = 

The energy eigenvalues are illustrated in Fig. 1. For Ml = lOto = H eV the circularly 
polarized light is in the off-resonance regime, where Aq = 0.6 eV [evA = 2.58 eV, top row 
of Fig. 1) or Aq = 0.8 eV {evA = 2.97 eV, bottom row of Fig. 1). The direct band 
gap of M 0 S 2 amounts to A = 0.83 eV and we have A = 37 meV, v = 0.5x10® m/s, and 
a = 3.193 A j^. For Aq = 0.6 eV the band gap is reduced to 0.46 eV for the iF'-valley 
(top right) and enlarged to 2.86 eV for the iF-valley (top left). The effect of the off-resonant 
light (Aq) can be tuned by varying the intensity, where various values have been achieved 
experimentally j^|. For Aq = 0.8 eV we observe that the spin splitting in the conduction 
band is increased and the band gap is reduced to 0.06 eV in the iF'-valley, so that only 
the iF'-valley is relevant, whereas for the TF-valley the band gap becomes 3.26 eV. In the 
following we restrict the discussion to the 77'-valley (p = —1). 


III. ORBITAL MAGNETIC MOMENT AND TEMPERATURE DEPENDENT 
ORBITAL MAGNETIZATION 


In order to study the Berry phase mediated thermoelectric transport, we consider the 


free energy, which, for a weak magnetic held B, is given by 

p *’-'=-1 ^ log (1 + ■ w --')) , 

P \ / 


( 6 ) 


where (3 = l/fc^T, k^ = 8.62 x 10“® eV/K is the Boltzmann constant, p is the Fermi energy, 
and T is the temperature. Eq. (6) can be simplihed by converting the summation into an 
integral, 

1 f cPk 


= - 


f3j (27ry 


(1 + ■ b) log (1 + , 


where 


= Vk X (Tp'(k) I^Vkl ^r'(k)) z 


^ U V k| 

is the Berry curvature. The energy Fp^(k) = — mp^(k) -B is modihed by the orbital 


(7) 


( 8 ) 


magnetic moment 


mr^(k) = 


-le 


n 


Vk^pAk) 

4 


H -E. 


s,-l 


Vk4/pAk))z 


(9) 








The orbital magnetization then is obtained as 


M*'^' = 



( 10 ) 


+ j-, I log (i+. 


where / is the Fermi distribution function. From Eqs. (4) and (5) we obtain for the z- 
component of the Berry curvature 

A — Aq + sX 


(k) = 

^ 2 [^ynky + (A - Ao + sA)2]3/2 

and correspondingly for the ; 2 -component of the orbital magnetic moment 


( 11 ) 


m 


s,-l 


(k) ^ -^El 


!iy(k). 


( 12 ) 


For hnite A — Aq + sX the orbital magnetic moment has a peak at A; = 0. For A = 0 
we obtain for A — Aq = 30 meV for a single valley and Aq = 0.8 eV an orbital magnetic 
moment of 35 Bohr magnetons. High magnetic moments have been predicted for systems 
involving orbital degrees of freedom, such as graphene 151. 

Using Eqs. (11) and (12) in Eq. (10), we obtain for T —)■ 0 for the conduction band the 
z-component 


^ 

2h 


1 - 


A — Aq + sA 
/i “h sA 


(13) 


which again can be enhanced by reducing the band gap by off-resonant light. Eq. (13) 
reduces to a previous result for gapped graphene in Ref. [151] in the limit of A = 0 and 
Aq = 0. As an example, for fi = 0.2 eV we obtain an orbital magnetization (by dividing 
the results of Eq. (13) by a typical layer thickness of 0.6 nm) of 0.1 Tesla, which is well 
detectable in experiments. 

In order to evaluate temperature effects, we study the spin and valley orbital magneti¬ 


zations ^ 


^ and M’' = ^ using Eq. (8), in Fig. 2 as a 


/■+!,—1 




function of the Fermi energy for temperatures of T = 160 K (left) and T = 360 K (right). To 
achieve spin and valley orbital magnetizations, respectively, we require A > |A — Aq| and 
A < |A — Aq|. is growing faster in the conduction band than in the valence band and 
is much/slightly larger than in the valence/conduction band. Increasing the intensity 
of the off-resonant light increases/decreases M'" in the conduction/valence band due to a 
slow/fast reduction of the band gap, which is due to the energy shift by SOC (see Eq. (4)). 
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As compared to Fig. 2, without off-resonant light we obtain about half the value for 
while is 100 times smaller (with opposite sign), since it is dominated by the band gap 
(the system is pinned to the valley transport regime). For T = 160 K we observe spin 
effects close to the Dirac point, whereas they are suppressed for T = 360 K. These effects 
will become clear in the next section, since the orbital magnetization is proportional to 
the Hall conductivity [l5|. Susceptibility measurements, electron paramagnetic resonance, 
x-ray magnetic circular dichroism, and neutron diffraction can be used to probe the orbital 
magnetization 28N30|. 


IV. SPIN/VALLEY THERMAL AND NERNST CONDUCTIVITIES 


Eq. (10) contains conventional (hrst term) and B erry phase mediated (second term) 
contributions. It has been demonstrated in Refs. IJ, ll5| that the conventional part does 
not contribute to the transport, while the Berry term directly modihes the intrinsic Hall 
current (which is obtained by integrating the Berry curvature over the two-dimensional 
Brillouin zone). In contrast to the Hall conductivity, the Nernst conductivity is determined 


not on 


y by the Berry curvature but also by entropy generation around the Fermi surface 


201. l21| and therefore is sensitive to changes of the Fermi energy and temperature. For 
a weak electric hel d E, the Hall current is given by jx = V^T) and the Nernst 


conductivity by 


cr*’ ^ = 

^xy 


eks 

h 


^ / ^«r‘(k)sr‘w, 

s,C ' 


(14) 


sylk) = -/(sy‘(k)-l<) ln/(£y‘(k)-^)-(l - /(£y‘(k) - ^)) in (1 - /(sy‘(k) - 1 *)) 
bein g th e entropy density. Recent experiments found that Eq. (14) describes graphene very 


well 


3lj. We obtain the transverse thermal conductivity 


_ kse ^ f (Pk 


D^’-^(k) <; - + - /i) (15) 


xy 


f)h (27r)2-< 


TT 


ris, —1 


5 ,- 1 / 


liT 


(l + - 2Li2 (1 - f{El’-\\E) - /i 


where Li 2 is the polylogarithm. Eqs. (14) and (15) can be simplihed in the limit of low 
temperature, using Mott relations 1^, to 


a. 


_ 


TT^ „ da: 


xy 


_ Jd rjn 


S, —1 

xy 




= —e 


dl^xy ^ 
dfi 


( 16 ) 
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and 


K 


s, — 1 
xy 


TT^ ^-2 


s,-l 
xy ) 


(17) 


where 

p 2 r A2 u 

= (U(E) - f.{E)) nr\k) (18) 

with f± representing the distribntion fnnction of the electron/hole bands. According to Eq. 
(16), the Nernst condnctivity is proportional to the derivative of the thermal condnctivity. 
We solve Eq. (18) for T = 0 by performing the integral and obtain in the case that the Fermi 
level is in the conduction band 


a!:-' = 


A - An + sA 


'xy 


‘2h jj, -\- sX 


(19) 


Eqs. (14) and (19) show that the spin (a® = a 






a 


J’ ^) and valley (a^ = ^ + 


a^y’ ^) Nernst conductivities are enhanced under off-resonant light. Without off-resonant 


light the spin Nernst conductivity is negligible because of the vanishing spin Hall conductivity 
in the limit A S> A (the system is pinned to the valley Hall regime). Due to the large band 
gap, the valley Nernst conductivity is also small for to < A. An enhanced spin splitting 
and corresponding giant thermoelectric transport in both the conduction and valence bands 
is achieved by reducing the band gap (to the range of A). The results in Eqs. (14) and 
(19) guarantee for monolayers of M 0 S 2 and related group-VI dichalcogenides an electrically 
tunable band gap to tailor the spin and valley transport. 

Experiments indicate that the thermoelectric properties can be understood by Mott rela¬ 
tions, which agree with experimental data for low temperature 3ll-l34|. In these experiments 
the thermoelectric properties, in particular the Nernst conductivity, have been measured for 
gapless graphene in a transverse magnetic held. The Nernst effect discussed in our work 
exists even without external magnetic held, being solely driven by the ehective magnetic 
held due to the Berry curvature. Note that Eq. (14) is more general than Eq. (16), because 
it goes beyond the linear temperature dependence. 

It has been found experimentally that the dependence of the thermoelectric transport on 
the gate voltage (Fermi energy) can be tuned by controlling the band gap in monolayer 

33| and bilayer 3^ graphene. Being the electrical response to the thermodynamic 


JJ, 


31 


perturbation, a giant thermoelectric transport is achieved when the bands come close to the 
Dirac point. In Fig. 3 (top) we show numerical results for ck®^, by evaluating Eq. (14), as 
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function of the Fermi energy at T = 160 K (left) and T = 360 K (right) and vary the band 
gap by off-resonant light as Aq = 0.8 eV (blue), Aq = 0.81 eV (green), and Aq = 0.82 
eV (red). For T = 160 K we hnd two peaks with negative values, where the lower peak 
is a consequence of the spin splitting and is washed out for T = 360 K. Our results show 
that a signihcant spin dependent transport can be observed in M 0 S 2 and related group-VI 
dichalcogenides at room temperature (or even above). Figure 3 (bottom) shows a^y as a 
function of the Fermi energy for T = 160 K (left) and T = 360 K (right) with = 0.6 
eV (blue). An = 0.65 eV (green), and An = 0.7 eV (red). We observe shifts of the peaks 
towards the Dirac point for increasing An, which reflects the reduction of the band gap. The 
amplitude grows for decreasing band gap. For An = 0.8 to 0.82 eV al.y is enhanced, whereas 
below An = 0.793 eV we are in the valley transport regime and obtain an enhancement of 
al.y. The transport is huge as compared to the case without off-resonant light. 

In general, it depends on the sign of the Berry curvature (compare Fig. 1) whether the 
Nernst conductivity is positive or negative. Our results are valid for elevated temperature 
in the experimentally relevant range 3l|. Furthermore, 7 ^ 0 when the Fermi energy 


is in the band gap, whereas Eq. (16) yields = 0 (being the derivative of which 

is quantized and independent of /i in this case). The demonstrated enhancement of the 
spin/valley transport due to the tunability of the band gap in M 0 S 2 and other group- 
VI dichalcogenides has been desired for thermoelectric applications since the discovery of 
graphene. Band gap opening by off-resonant light has been achieved in Ref. 23| for the 


surface states of topological insulators and can also be used for monolayer M 0 S 2 , since 


transistors [5j and amplihers already have been realized 


V. CONCLUSION 

We have derived analytical results for the thermoelectric transport in monolayer M 0 S 2 
and related group-VI dichalcogenides in the presence of off-resonant light. We have shown 
that an increased intensity of the light reduces the direct band gap and results in a strong 
spin splitting in the conduction band and, therefore, in dramatic enhancement of the ther¬ 
moelectric transport. The spin splitting in the conduction band (in contrast to the valence 
band) is negligible without external perturbation (such as off-resonant light). The band 
gap even can be tuned to zero with giant spin splitting in both the valence and conduction 







bands. The enhancement of the spin/valley transport properties demonstrated here is de¬ 
sired for spin/valley dependent thermoelectric devices, whereas the tunable band gap opens 
new directions for fundamental transport experiments. 


Appendix A 

The time dependence in Eq. (1) can be understood as the sum of two second-order virtual 
photon processes, where electrons absorb and then emit a photon and electrons emit and 


then absorb a photon 


22| . Within Floquet theory, we use the fact that 


=-\og[U{To)] 


with 


U{To) = texp 


fTo 

-i / 

Jo 


^ exp[-*i7-;To], 


where T is the time ordering operator. We consider the Fourier decomposition 


(A.l) 


(A.2) 


—iflt 


n=—oo 

where we have used 

1 fTo 

^0 Jo 

Expanding the exponential function in Eq. (A.2) in a Taylor series yields 


exp 


fTo 

-i / 

Jo 


(A.3) 


(A.4) 


fTo I—A 2 r'^o 

= l-i + H^^’^{t2)dt2 + ... 

Jo 2 Jq Jq 

(A.5) 


and thus 


17(ro) ~ t (^1 - ^ ^ ° ^ ° H^^'\ti)dti ^ ° W^’\t2)dt2 ) . (A.6) 

By applying the time ordering operator in Eq. (A.6) we arrive at 


t/(To) ~ 1 - i / ° H^^'^{t)dt -I f ° dti [ ' dt2H^^'^{ti)H^^’'^{t2) (A.7) 


./o 

1 


0 T 2 


I dt2 / dtiH^^’\t2)H^^’^{ti). 

0 Jo 
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Executing the integration in Eq. (A.7) with the help of Eq. (A.3) and reordering the terms, 
we obtain 


17(To) ~ 1 - ^ H^.f] 


(A.8) 


= I - --{H:y^)X 


Tn" + - 


and thus 


Kff = nt” + ^ Kf] - im-”. n‘Jf]] + X*"'’ 

with = vijjaxPx + (^yVy) + Aa^ — XrjSzCTz + XrjSzi which describes a static honeycomb 

lattice with hopping to (in a standard tight binding notation) and a band gap of 2A. When 
Ml ^ to and A then the off-resonance condition is satished and perturbation theory can be 
applied. Since the first order term in Eq. (A.9) vanishes, we simplify the second order term 
to arrive at the effective Hamiltonian of Eq. (3). 
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FIG. 1. Band structure of monolayer M 0 S 2 in the presence of off-resonant light and SOC for 
the iL-valley (left) and iF'-valley (right) using HQ = lOto = 11 eV, A = 0.83 eV, A = 37 meV, 
V = 0.5x10^ m/s, and a = 3.193 A. The top row refers to Aq = 0.6 eV and the bottom row to 
An = 0.8 eV. 
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FIG. 2. Spin (top) and valley (bottom) orbital magnetizations (divided by the layer thickness 0.6 
nm) as a function of the Fermi energy for T = 160 K (left) and T = 360 K (right). For we 
use Afi = 0.8 eV (blue), Aq = 0.81 eV (green), and Aq = 0.82 eV (red), whereas for M’' we use 
An = 0.6 eV (blue), An = 0.65 eV (green), and An = 0.7 eV (red). All other parameters are the 
same as in Fig. 1. 
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FIG. 3. Spin (top) and valley (bottom) Nernst conductivities as a function of the Fermi energy 
for T = 160 K (left) and T = 360 K (right). For a%y we use An = 0.8 eV (blue), An = 0.81 
eV (green), and An = 0.82 eV (red), whereas for we use An = 0.6 eV (blue). An = 0.65 eV 
(green), and An = 0.7 eV (red). All other parameters are the same as in Fig. 1. 
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